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With M(t) := supj,g[g A{s) — s denoting the running maximum 
of a fractional Brownian motion A{-) with negative drift, this paper 
studies the rate of convergence of P(M(i) > x) to P(M > x). We define 
two metrics that measure the distance between the (complementary) 
distribution functions f{M{t) > ■) and P(Af > ■). Our main result 
states that both metrics roughly decay as exp(— where i9 is 
the decay rate corresponding to the tail distribution of the busy pe- 
riod in an fBm-driven queue, which was computed recently [Stochastic 
Process. Appl. (2006) 116 1269-1293]. The proofs extensively rely on 
application of the well-known large deviations theorem for Gaussian 
processes. We also show that the identified relation between the de- 
cay of the convergence metrics and busy-period asymptotics holds in 
other settings as well, most notably when Gartner-Ellis-type condi- 
tions are fulfilled. 



1. Introduction. Let (A(s))sgM be a centered fractional Brownian mo- 
tion (fBm), that is, a stochastic process, such that for aU s G R, A[s) obeys 
a Gaussian distribution with mean and variance |sp^, for H G (0, 1). fBm 
has recently become one of the key models in the applied probability litera- 
ture, because of a number of interesting features. For H ^ (^) 1)) correlations 
decay so slowly that the process qualifies as long-range dependent; choos- 
ing H € leads to positive correlations, whereas H € (0,^) results in 
negative correlations. Also, fBm exhibits selfsimilar behavior, in that A{as) 
has the same distribution as A{s). Its use has been advocated in several 
practical settings; see, for example, [21, 28] for networking applications; for 
applications in finance, see, for example, [2]. 
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Motivated by these applications, substantial attention was paid to the 
analysis of regulated fBm, or fractional Brownian storage [27]. With the 
storage process defined through 

Q{t) := sup A{t) - A{s) -{t- s), 

s<t 

distributional properties of the steady-state storage level Q := lim(_>oo Qit) 
can be used to describe the performance of a network element. The sta- 
tionary storage level Q being distributed as M := supg^QA{s) — s, see for 
instance [22], Sections 5.1 and 5.2, a considerable amount of work [12, 16, 
25, 26] has been devoted to the characterization of the distribution of M, 
that is, the supremum attained by an fBm with negative drift. The results 
obtained are predominantly asymptotic in nature; most notably, an explicit 
function ip{-) was identified [16] such that 

P(M > x)/ip{x) 1 

as X — > oo. Clearly, ip{x) can serve as an approximation of ¥{M > x) for 
large x; we lack, however, accurate approximations or bounds for small or 
moderate values of x. 

As mentioned above, asymptotic results are available that describe P(Af > 
x) for X large, but considerably less is known about the convergence of the 
running supremum 

M{t) := sup A{s) - s 

se[o,t] 

to its limiting distribution M. The primary goal of the present paper is to 
determine the speed of this convergence. Knowledge of this speed of conver- 
gence is useful in several contexts. In the first place it provides information 
on the question whether at the time scale at hand, that is, t, it is justified 
to approximate P(M(t) > x) by P(M > x). Loosely speaking, if the conver- 
gence to the limiting distribution is fast (i.e., a short relaxation time), then 
such a procedure gives an accurate approximation, whereas in case of slow 
convergence transient results need to be used. Supposing that one wishes 
to estimate P(M > x) [and hence also P{Q > x)] by performing a stochastic 
simulation, the speed of convergence also provides us with useful insights 
with respect to the design of this simulation experiment. It, for instance, 
sheds light on the question of whether it is more efficient (in terms of the 
variance of the resulting estimator) to have one long run or multiple shorter 
runs. In addition, it yields guidelines as to how long one should simulate 
A{s) — s in order to be able to accurately estimate P(M > x), that is, to find 
an appropriate simulation horizon t such that P(M(t) > x) approximates 
P(M > x) sufficiently well (in terms of some explicitly specified criterion); 
cf. [11]. 
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In principle, the rate of convergence (as t grows large) of P(M(t) > x) to 
P(M > x) could depend on the value of x. Obviously, one could use several 
distance measures, each incorporating this dependence on x in a specific 
way. First notice that 

P(M >x)- P(M(t) >x)= P(M > x, M{t) < x) =: 7(3;, t) > 0, 

as the event {M{t) > x} implies {M > x}. Two possible distances are the 
following: 

• Kolniogorov-Smirnov distance (sup-norm). We define (for random vari- 
ables X and Y) di{X,Y) by sup^ \F{X < x)-F{Y < x)|. It is weh known 
that di is a distance. The first metric used in our paper is based on the 
distance di: 



Bi(t) measures the maximum distance between the distribution functions 
of M and M{t). 

• Integral distance (Li-norm). It is a well-known fact that d2{X,Y) := 
\¥{X < x) — ¥{Y < x)\dx is a distance, too. The second metric con- 
sidered in our paper is based on the distance d2- 



D2(i) measures the total distance between the distribution functions. Note 
that ©2(0 can be interpreted as EM - EM{t). 

The goal of the paper is to identify the asymptotics of the Bi(t) for t large, 
i = 1,2. Our main result is that the decay rates of both distance measures 
coincide, and are equal to asymptotics of the busy-period distribution in an 
fBm-driven queue, which were recently identified in [23]. 

The structure of this paper is as follows. In Section 2 we recapitulate a 
series of results on the large deviations of fBm, most notably (the generalized 
version of) Schilder's theorem. We also recall the main results on busy- 
period asymptotics [23], which enable us to state the main results of our 
paper. Section 3 presents a number of auxiliary results that are used in 
Section 4 (in order to determine the asymptotics of the Kolmogorov-Smirnov 
distance) and Section 5 (in order to determine the asymptotics of the integral 
distance). In Section 6 we consider the situation of short-range dependent 
input (or, more precisely, the situation in which so-called Gartner-Ellis-type 
conditions are met), to show that also in this regime the asymptotics of the 
Oi{t) are equal to those of the busy-period distribution. Section 7 concludes; 
it includes a procedure for determining the simulation horizon. 



Bi(t) := di (M,M(t)) = sup7(x,f). 



x>0 
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2. Preliminaries and main results. In this section we recall a number of 
useful results from the literature. Emphasis is on busy-period asymptotics 
recently identified in [23]. We then state our main results. 

2.1. Generalized Schilder. Informally, the generalized version of Schilder's 
theorem provides us with a "rate functional" I(-) such that 

Pn[y] :=P(^^ G PS exp(-n jnf I(/) 

In other words: in this large-deviations setting, the probability of interest 
decays exponentially in n. The "~" in the above statement should be inter- 
preted as follows: under mild conditions on the set ,y (more concretely, if 
^ is an I-continuity set), the decay rate oi pnl-^] is given by 

lim -logpr,[y]=- inf 

n^oo n fey 

Apart from the Brownian case H = ^, the "rate functional" cannot 
be given explicitly. It is defined through 



^■^■^ loo, otherwise, 

where ^ is the reproducing kernel Hilbert space related to the process A{-) — 
see for details [1, 9]. Here ||/|| := ^/ (/, /), where {f,g) is a suitably defined 
inner product between f,g It is noted that can be interpreted as 
a measure for the "likelihood" of a path /; the path / = is the sole path 
that gives I(/) = 0, while for other paths / the rate is strictly positive. 
The path /* := arginfjg,j^I(/) is usually called the most likely path in large 
deviations literature, and it has the interpretation that, conditional on the 
fBm being in the set =y, with overwhelming probability it will be close to 

/*; cf. [10]. 

Later in this paper we repeatedly use the following property. Suppose 
/ E and g is defined by g{r) = af{f3r), for a, /3 > 0. Then 

(1) ||5||=a/3^||/||. 

For the purposes of the present paper, more background on "generalized 
Schilder" is not required; see for a complete account [1, 9, 23]. 

2.2. Busy-period asymptotics. In [23], and its predecessor [29], the focus 
was on computing the asymptotics, for large values of t, of P(A' > t), where 

K := inf{t > : Q{t) = 0} - sup{t < : Q{t) = 0} 
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is the ongoing busy period at time 0. In [29] it was shown that F{K > t) de- 
cays roughly in a Weibulhan way, that is, as exp(— for some positive 
constant "i?. More precisely, it obeys the following logarithmic asymptotics: 

(2) lim -rr^ logPfi^ >t) = -^ where •& := inf 1(f). 

i— >oo f££lg 

Here ^ is the set of paths that remain above the diagonal on the interval 
[0,1]: 

^ := {/ G ^ I Vr e [0, 1] : /(r) > r}. 

Note that the set ^ can be regarded as the set of feasible paths that 
correspond to an intersection of events (reflected by the V-quantor). Where 
unions are usually easy to deal with, finding the minimizing path in an 
intersection is typically hard (although often rather precise bounds can be 
found; see for instance [24]). In [23] we succeeded in determining the right- 
hand side of (2), as well as the corresponding minimizing path /* £ 

In this particular setting, the most likely paths turn out to have a re- 
markable shape. For H > ^, the most likely path is at the diagonal in some 
interval [0, s*], and also at time 1, but strictly above the diagonal in between. 
For H < ^, the corresponding path departs immediately after time from 
the diagonal, but returns to it strictly before time 1 and continues along it 
until time 1 — see Figure 1. The corresponding decay rate is given in [23], 
Theorem 24. 

2.3. Main results. We now present the main results of this paper. They 
entail that both ID'i(t) and ID)2(i) decay as the probability ¥{K > t) of the 
busy period exceeding t. 

Theorem 2.1. (i) For the Kolmogorov-Smirnov distance we have 

}}^l^^o,B,it) = -in^Iif) = -^. 
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(ii) For the integral distance we have 

}}^ ^ = - inO(/) = 

Part (i) and (ii) of this theorem will be proven in Sections 4 and 5, re- 
spectively, whereas Section 3 develops a number of useful tools. 

3. Auxiliary results. In this section we derive a number of results that 
are needed to prove Theorem 2.1. The following alternative expression for 
7(x,i) turns out to be useful. 

Lemma 3.1. For any a; > 0, t > 0, 

7(x, t) = P(Vr e [0, t] : A{r) < x + r; 3s > t : A{s) > x + s) 

VrG[0,l]:^<|+r;3s>l:^>|+. 

Proof. The first equality is a matter of rewriting {M > x, M{t) < x} in 
terms of the process A{-). The second equality follows from the self-similarity. 
□ 



In the sequel we extensively use the following sequence of probabilities: 
p„(5):=p(^VrE[0,l]: 



^^'^hs + r;3s>l:m>S + s 



in \jn 
for 5 > 0. We also define their exponential decay rates by 

J{S) := lim -logp„(5), 

again for > 0; here J(0) denotes the limit of J(S) for (5 | 0. 
Define as the paths / in the set 

(3) i4 := {/ G ^ I Vr G [0, 1] : /(r) < 5 + r; 3s > 1 : /(s) > 5 + s}; 

also := s^Q. The proof of the following result is an immediate consequence 
of (the generalized version of) Schilder's theorem [1, 9]. 

Lemma 3.2. For any 6 >0, 

(4) J(<^)=- inf I(/). 
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=c/ := {/ G ^ I Vr G [0, 1] : /(r) < r; /(I) = 1}; 
^ := {/ G ^ I Vr G [0, 1] : /(r) > r}. 

The following result concerns a translation of J(0) in terms of our previous 
result on busy periods, as mentioned in Section 2.2. 

Proposition 3.3. 

infJ(/)= infJ(/)= inf 



Proof. Due to continuity arguments (cf. the proofs in Section 4 of [29]) 
the decay rate corresponding to the most likely path in £^ is the same as 
that of the most likely path in 

{/ G ^ I Vr G [0, 1] : /(r) < r; 3s > 1 : f{s) > s}. 

Let s* be the smallest s > 1 such that f{s) > s. Define the path / through 
/(r) := f{rs*)/s*. Then /(I) = 1 and, due to (1), 

||/||2 = (,^)2//-2||/||2<||J||2_ 

This implies the first equality. 

A time shift argument trivially gives that inf is equal to inf 

with 

=£/_:={/ G ^ I Vr G [-1, 0] : /(r) < r; /(-I) = -1}. 

Now reverse time, and we obtain that this infimum is also equal to inf 
with 

J := {/ G I Vr G [0, 1] : /(r) > r; /(I) = 1}. 

The analysis of [23] implies that the infima over ^ and ^ coincide, which 
proves the stated. □ 

To obtain a better intuitive understanding of Proposition 3.2, the reader 
may compare the most likely busy-period path (i.e., the "cheapest path" 
in as depicted in Figure 1, with the most likely path in = s^q, as 
depicted in Figure 2. 

Lemma 3.4. (i) With, for 6>0, 

(5) ^5:={/G^|VrG [0,l]:/(r)>r;/(l) = l + 5}, 

infjgg^^ II/IP increases in 6, for 6 G [0,H~^ — 1]. 
(ii) With, for 0<e<l and 5>0, 

&5,e-={f^^\yre[0,l]:f{r)>r-e;f{l) = l + S-e}, 

infjgg^^^ II/IP increases in 5, for 6 £ [0,H~^ — 1]. 
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Proof. In order to settle claim (i), first observe that the probability 
qn{6) := pfvr G [0, 1] : ^ > r; 3s > 1 : ^ >5 + s 

\ x/lT' \/TI 



decreases in 6, so consequently also 
(6) lim - logqn{5) = -I- inf 

decreases in 6; the equality is due to (the generalized version of) Schilder's 
theorem, and 

#5 := {/ G ^ I Vr G [0, 1] : /(r) > r; 3s > 1 : /(s) > <5 + s}. 

Pick an arbitrary path / in this set, and let s* be the smallest s > 1 such 
that /(s) > s + Then consider the path /, defined by 

/»-(^) /("*)■ 

Note that /(I) = 1 + 6, and, because / lies in the set ^s, ^or all r £ [0, 1], 

it is easily verified that the right-hand side of the previous display is at least 
r. We conclude that / is in &s as well. Moreover, as before, any path in 
can be replaced by a path in with a smaller norm: due to (1), 

ii/T=(^)'(^^ni/f <ii/f. 

Here it is used that for 6 G [0, - 1] and s > 1 it holds that (1 + 5)s^ < 
s + 6. In other words, for these 6 we can replace the set &s in (6) by 
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Hence 

inf 

decreases in S as well, which implies claim (i). 
The proof of claim (ii) is similar. The probability 

qn{6, e) := P fvr G [0, 1] : ^ >r-e;3s>l:^>6 + s-e), 

decreases in 6. Again, for any / in 

^5,e-={f^^\^re[0,l]:f{r)>r-e;3s>l:f{s)>6 + s-e}, 
and s* as defined before, we can define 



which has a smaller norm than / for any 5 G [0,H~^ — 1], and which lies 
in &s,£] the latter statement follows from f{l) = l + 5 — e, in combination 
with, for r G [0, 1] , 

where the first inequality is due to / G ^s,e i ^-nd the second inequality due to 
the fact that, obviously, {s* — l){Sr + e(l — r)) > (realize that 1 + 5 — e> 
because e < 1). We have now established claim (ii). □ 

We now establish a useful lemma on the behavior of J (6), that is, the 
decay rate of Pn{5), as a function of 5. 

Lemma 3.5. (i) J{5) decreases for 6 G {0,H^^ — 1]. 
(ii) J{5) < J(0) for 6 > H'^ - 1. 

Proof. First consider part (i). Due to continuity arguments, and "gen- 
eralized Schilder," the decay rate J[5) equals the decay rate corresponding 
to the most likely path in 

(7) {/ G ^ I Vr G [0, 1] : /(r) < r + 5; 3s > 1 : f{s) >s + 5]. 

With arguments identical to those used in the proof of Lemma 3.4, we can 
show that for 5 G [0, — 1] we can replace the set (7) by 

(8) {/G^IVrG [0,l]:/(r)<r + 5;/(l) = l + 5}. 

Reversing time, we observe that the decay rate corresponding with the 
most likely path in (8) equals the one corresponding to the set as defined 
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through (5). Observe that the event corresponding to the latter set becomes 
increasingly rare when 6 grows; Lemma 3.4(i) now implies claim (i). 
Now consider the second claim. Clearly 

J(5) < lim - logpf 3s >1: >5 + s]. 

The decay rate in the right-hand side equals (see, e.g., [22], Exercise 6.1.3) 

which reduces, for 5 > H~^^ — 1, to 

(9) 1/ ^ /I 



2\1-HJ \H 

this is a decreasing function, with value —-^H'^ for 6 = — 1. Conclude 
that, for 5 > - 1, J{6) < It is left to prove that J(0) > -\H-^ . 

To this end, we use that in [29] it was shown that 

-\-AH) < J(0)<-i; 

.jr. 1 r(3/2-g) 

^' H{2H -1){2-2H) T{H -1/2)T{2-2H)' 

As seen from Figure 4.1 in [29], we need to check <f{H) < i/^'^ for H in 
the neighborhood of 1. Observe that both functions have value 1 for H = 1. 
As ip{-) is concave in this neighborhood, and convex, we need to verify 
whether ip'{l) > —2. Calculations yield that, with 7em the Euler-Mascheroni 
constant 0.5772, 

r'(i/2) 

<^'(l) = -3 - 2 • - 27EM « -0.22 > -2. 

This proves claim (ii). □ 

Lemma 3.5 is already a first indication that Theorem 2.1(i) indeed holds, 
as seen as follows. Observe that (provided that the limits exist), using 
Lemma 3.1, with a defined as (2 — 2H)~^, 

= lim 2-2H log sup 7(^1^) 

1 / Air) X 

= lim -logsupP Vr G [0,1]:— ^ < — + r; 
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1 / A(r) 

= lim -logsupP Vr G [0,11:—^ <x + r; 

n^oon x>0 \ 

3s >l:^pJ- >x + s] 
Vn J 

= lim — log sup pn(x). 

n^oo n x>0 

In other words, if one can interchange the hmit and supremum, then 
J™ ^2 iQglPi(^) = sup J(a;), 

which equals J(0) due to Lemma 3.5; applying Lemma 3.2 and Proposition 
3.3, conclude that this would also mean that Theorem 2.1(i) holds. The 
goal of Section 4 is to prove that the limit and supremum can indeed be 
interchanged. 

In light of the fact that, in a large-deviations setting, the decay rate of an 
integral is, under rather general conditions, determined by the decay rate of 
the maximum of the integrand, it is now also expected that indeed Theorem 
2.1(ii) holds. Section 5 is devoted to substantiating this claim. 



4. Proof for Kolmogorov Smirnov distance. 



Proof of Theorem 2.1(i). As seen above in (10), 

lim ^ logBi(t) = lim - logsupp„(a;), 

implying that the lower bound is trivial, as for all e > 0, 

^lim ^2iH^logDi(t)> J(e). 

Now let e i 0, and apply Lemma 3.2 and Proposition 3.3. The upper bound 
is proven in the following steps: 

(i) First observe that, taking for convenience M as a multiple of e > 0, 



/M/e \ 

supp„(x) < y^Pnik) + supp„(a 

x>0 Vfc=l / =^>M 



(11) 

withp^(A;) :=sup^g((fc_i)£_fc£]Pn(2;) 
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(ii) Trivially, 

sup pn{x) < sup f(3s>1: > x + s] =f(3s > 1: >M + s 
x>M x>M \ \Jn J \ 

Also, as seen in (9), for M > — 1, 



n^co n V \/n J 2\1 — HJ \H 



We saw in the proof of Lemma 3.5 that, for M > H ^ — 1, this expression 
is majorized by J(0), and hence 

(12) lim - log sup Pn{x) !^ J (0). 

n x>M 

From now on we pick M : = — 1. 
(iii) Also, 

plik) < pfvr G [0, 1] : < fee + r; 3s > 1 : > {k - l)e + s] ; 

V \/n Jn J 



because of "generalized Schilder," we thus obtain 

lim ilog/,(A:)<- inf I(/), 

where 

:= {/ G ^ I Vr G [0, 1] : /(r) < x + r; 3s > 1 : /(s) >x-e + s]. 

Applying Lemma A.l in the Appendix to the right-hand side of (11), and 
using (12), we obtain 

(13) lim i logsupn„(a;) < max/ max (— inf I(f)Vj(0)]>. 

(iv) We now show that, for x = e,2e, . . . , M, 

(14) - iiif !(/)<- irif I(/). 

This is done as in Lemma 3.5(i). First, using the arguments of the proof of 
Lemma 3.4, we can restrict ourselves for G {e, 2e, . . . , M} to the paths that 
attain the value 1 + x — e at time 1, that is, paths in 

{/ G ^ I Vr G [0, 1] : /(r) <x + r; /(I) = 1 + x - e}. 

Reversing time yields that this is equivalent to finding the most likely path 
in 

{/ G ^ I Vr G [0, 1] : /(r) > r - e; /(I) = l + x-e}. 

We now observe that the event corresponding to this set is increasingly rare 
for growing x, as made precise by Lemma 3.4(ii) (where we have chosen 
e < 1). Hence, for x = e, 2e, . . . , M, claim (14) now follows. 
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(v) Hence the right-hand side of (13) is bounded from above by 
(15) ™ax|-/^f I(/),'/(0)|. 

By letting e j 0, 

MO) := (- inf I(/)) i (- inf I(/)) = J(0). 

Now the stated follows from letting e | in (15), and application of Lemma 
3.2 and Proposition 3.3. □ 

5. Proof for integral distance. 

Proof of Theorem 2.1(ii). We start by estabhshing the lower bound. 
Evidently, for s > arbitrarily chosen, and, as before, q := (2 — 2H)~^, 

/ j{x,t)dx> / j{x,t) dx > en°' ( inf j{x,t)). 

Jx>0 J x&[en°- ,2en°'] Vx6[en",2en"] / 

Hence, for e > arbitrarily small, using ■ logn — > 0, 

1 1 

lim inf logB2(t) > liminf — log( inf j{x,n" 

t^oo t^oo n ^xe[en°' ,2en°'] 



liminf — logf inf pn(x)\ 



n ^a:e[e,2e] 

It is straightforward that 

inf pn{x) > pfvr G [0, 1] : < e + r; 3s > 1 : >2e + s 

x&[e,2e] V yri 



The lower bound now follows from letting e | 0, together with the usual 
continuity arguments and time reversal. 

We now turn to the upper bound. Obviously, for any e > 0, 

» oo 

/ ^{x,n°')dx <n°'^'^i sup 7(x,n")) + V" ( sup 7(x,?i") 

"'a:>0 ^xe(0,n-+=] ^ fc= J ^£[^''^+1) 

We consider the decay rates (in n) of both terms. First focus on the first 
term; because of, again, ■ logn 0, and Theorem 2.1(i), 

lim sup — log (n°'^^ ( sup 7 (x , n° 



n— >oo n 

(16) 



< lim sup — logf sup 7(x,n")) =— inf I(/). 

n^oo n ^x6(0,oo) ^ /e^ 
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Also, applying part (i) of Lemma A. 2 in the Appendix, there are positive 
constants k, A such that 



1 / °° ^ 
limsup — log ( sup 7(x,n") 

n— >oo n \,i '' 



xG[fc,fc+l) 



(17) 



< limsup -log ^ ( sup P(M>x))) 

'^^^ " Vfc=Ln-+-J^6[fc,fc+l) V 

< limsup -log ^ ( — ""^ \^2-2/f^ 



sup Kexp(— Ax 

j7^+ej ^X6[fc,fc+1) 
-j^ / oo \ 

<limsup — log ^ Kexp(— AA;^~^^) . 



n— >oo ?T- 



VA;=[n"+Ej 



Now part (ii) of Lemma A. 2 can be applied: there exist positive constants 
K,X such that (17) is bounded above by 

limsup - log(Sexp(-A(n"+")2-2^)) = -oo, 

n— »oo TL 

where the last equality follows by recalling that q = (2 — 2H)~^ . Combining 
this with (16) and using Lemma A.l completes the upper bound. □ 

One might think that asymptotic results for Di are stronger than those for 
5 as they involve the whole distribution rather than just the first moment. 
However, the metrics are different and a priori the asymptotics could differ. 
Consider for instance a situation in which P(Z > x) = x~" (where x > 1, 
a > 1) and F{Z{t) > x) = l(i^()X~" (again x > 1, a > 1), where we have that 
di{Z, Z{t)) = t"" but d2{Z, Z{t)) = t-"+V(a - 1). It is noted however that, 
despite the fact that the asymptotics may differ, both notions are strongly 
related under rather broad circumstances; see, for instance, the proof of 
Theorem 2.1(ii), where we could make use of the result for the asymptotics 
of Di(t) to establish the corresponding result for the asymptotics of D2{t). 

6. Analogous results for short-range dependent input. The analysis for 
fBm shows that the logarithmic asymptotics of both distances coincide with 
those of long busy periods. One may wonder whether such a property is valid 
under more general circumstances. One could pursue to extend the class of 
models for which this result holds to Gaussian processes with regularly vary- 
ing variance functions; cf. [8, 10]. In this section we focus on non-Gaussian 
processes, namely, short-range dependence processes that obey Gdrtner- 
EUis-type conditions; see, for example, [17]. 
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To this end, with A{t), as before, the traffic generated by a process with 
stationary increments in a window of length t (which we can assume to have 
zero mean, without loss of generality), consider for x > 

K{x) := lim - logP(Vr G [0, t] : A{r) < xt + r;3s > t: A{s) >xt + s). 

t— >oo t 

It is trivial to rewrite this decay rate to 

lim - logP(Vr G [0, 1] : • A{rt) < x + r; 3s > 1 : • A{st) >x + s). 

i— »oo t 

First define the (asymptotic) cumulant function 

A(s) := lim -logEe'^^, 

t^co t 

which we assume to exist; this essentially means that the input traffic is 
short-range dependent. Then one can define the large deviations rate func- 
tion by its Legendre transform 

I{a) := sup(sa — A(s)). 

s 

It is readily verified that, under mild conditions, the decay rate in the pre- 
vious display is bounded from above by 

(18) lim - logF{3s>l:t~'^ ■A{st)>x + s) = -inf si 

t— >oo t S>1 

these mild conditions in particular relate to the behavior of the input process 
between grid points, as formalized in [12], Hypothesis 2.3. 

We now study under which conditions this upper bound is actually tight. 
First observe that for any T > 1 , 

K{x) > lim - log P(Vr G [0, l]:t~^ ■ A{rt) <x + r- 

t^QO t 

3s G (1, T] : • A{st) >x + s); 

below we specify how T should be chosen. Let the infimum in the right-hand 
side of (18) be attained in > 1. Suppose t~^ ■ A{4) satisfies a sample-path 
large deviations principle (sp-LDP) on [0,T] of the Mogulskii type, with rate 
function I(-) and supremum-norm, then 

K{x) > -!(/) 

for any feasible path /, that is, all / in as defined by (3); here !(/) := 
/q I{f'{t)) dt. Now verify that the path given by /*(s) = s ■ {x + s*) / s* for 
s G [0, s*], and /(s) = x + s* for s > s*, is indeed feasible; also 

i(r)=s^/'''^+'' 
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SO that we can conclude that indeed 



Kix) = -mfsl( 

s>l \ S J 



Now that we have an expression for K{x), we wonder whether we can prove 
the decay of the Oi(t) is similar to that of the tail distribution of the busy 
period. To this end, first note that K(x) is a decreasing function of x. It is 
relatively straightforward to check that this fact entails that, for i = 1,2, 

lim -logDi(t) = K(0) = -1(1). 

This argumentation indicates that the decay rate of Oi(t) indeed coincides 
with the busy-period asymptotics —/(I), like in the fBm case, as long as 
an sp-LDP is available. For a broad class of discrete-time processes satis- 
fying a Gartner-Ellis condition (covering discrete-time Markov modulated 
processes), such an sp-LDP was proved by Chang [6], whereas for Levy pro- 
cesses see, for example, [7]. It is clear that, to make the above argumentation 
work, it is sufficient that T is chosen larger than s*. 



7. Discussion and concluding remarks. 



Decay of the correlation function. In Section 6 we showed for short-range 
dependent input that, under specific regularity conditions, Bj(t) (i = 1,2) 
decay essentially exponentially in t, and this decay roughly coincides with 
that of the tail of the busy-period distribution. Let us now consider the 
asymptotics of the covariance Cov {Q{0),Q{t)), with, as before, 

Q{t) := sup^(t) - A{s) -{t- s), 

s<t 

assuming that the queue is in equilibrium at time 0. In [14] it was shown 
for short-range dependent Levy input (with no negative jumps), that also 
Cov {Q{0),Q{t)) has the same asymptotic behavior as the tail of the busy- 
period distribution: for i = 1,2, 

lim -logCov(Q(0),Q(t)) = lim - logDi(t) = -/(I) 

t— ►OO t t—fOO t 

with /(I) as defined in Section 6. 

Our main result states that also in the fBm case, we saw that the busy- 
period asymptotics and those of Di{t) {i = 1, 2) match. In light of the above 
findings for short-range dependent Levy input, this suggests that also in the 
fBm case 

lim ^^logCov(Q(0),Q(t)) = 

t—>oo t 
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Based on the recent results in [13], however, we expect that this is not true. 
Instead, we anticipate that the asymptotics of Cov{Q{0),Q{t)) are roughly 
polynomially, or, more precisely, decaying as t^^~^ , which is equally fast as 
the asymptotics of Cov(y4(0, l),A{t,t + 1)). A formal proof of this property 
is still lacking, though. As we feel that the determination of the correlation 
asymptotics of reflected fBm is an important open problem, we state it as a 
conjecture. 

Conjecture 7.1. For some constant 7 e (0, c3o), 
lim t^-"^" CoY{Q{Q),Q{t))=-i. 

t— >oo 

Relation with busy-period asymptotics. The relation between the asymp- 
totics of the metrics considered and those of the tail of the busy period has 
been observed before in various contexts before. It is noted that in a number 
of specific Markovian settings one has shown that the rate of convergence 
does not depend on the initial state, which is a result that is reminiscent 
of ours; we refer to [4, 5, 30], which are closely related to generic results 
by Kingman [19, 20], for the case of birth-death processes. Regarding the 
fact that the decay rate of the Di(t) coincides with that of the tail of the 
busy-period distribution, we mention a classical result by Kingman [18], 
Theorem 7, Lemma 7, in the setting of an M/G/1 queue, and results for 
various queueing systems by Blanc and van Doorn [3]. For the special case 
of spectrally positive Levy input the double Laplace transform of 7(2;, t) can 
be determined in a similar way as in [14], and inversion techniques can be 
applied to identify the exact asymptotics [i.e., a function ijjx{-) is found such 
that il^x{t) ■ 7(3;, t) — > 1 as t — > 00]; elementary calculations show that then 
also the busy-period decay rate appears. 

On an intuitive level it is clear that there is a relation between the rate of 
convergence to stationarity and the tail of the busy-period distribution. In 
case of short-range dependent input (and in particular when busy periods 
are independent) the autocorrelation essentially breaks when the busy pe- 
riod ends. Our result shows that this heuristic carries over to the setting of 
long-range dependent input (in which there is still substantial dependence 
between consecutive busy periods). We remark however, that under gen- 
eral conditions (i.e., irrespective of the traffic being long-range dependent or 
short-range dependent) upper bounds on the metric Oi{t) in terms of the 
tail distribution of an ongoing busy period, can be constructed relying on a 
coupling argument [15]. 

Use of convergence estimates in simulation. If one aims at estimating the 
probability P(Af > x) through simulation, one needs to truncate the infinite- 
time horizon to some finite value, say t. This evidently always implies an 
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underestimation. Obviously, one needs to choose t sufficiently large such that 
the error made is negligible. Let Tx be defined as the smallest s such that 
A{s) - s = x; then it holds that P(M > x) = P(T^ < oo). Then one could, 
for instance, require for small e > (for instance 5%) that t be chosen large 
enough that 

P(t < Tx < oo) 
F{Tx < cx)) ^ ^' 

The numerator equals 7(x, t), and can therefore be approximated by exp(— ■!? x 
^2~2H\^ (noting that Lemma 3.5 suggests that this is a conservative estimate), 
whereas the denominator can be bounded from below by [27] 1 — ^>(x + 1*), 
with $(•) denoting the standard normal distribution function, and t* : = 
xH/{l — H); 1 — ^{x + t*) is approximated by 



In this way we can compute an estimate for the simulation horizon t: 



t> 



loge 1 f X \2-2^^/ 1 l/(2-2H) 



i9 2^\l-H \H 



The horizon grows in x, as expected, and it does so in a linear fashion for 
X large. A procedure for the discrete-time counterpart is detailed in [11], 
see also Proposition 8.1.1 in [22]; it shows the same qualitative behavior as 
a function of x. A procedure for numerically computing •& can be found in 
[23]. 

As mentioned in the Introduction, a second evident application of our 
estimates relates to another issue in the design of the above simulation 
experiment: it provides insight into the question whether it is more efficient 
to simulate one long run, or to simulate multiple shorter runs. 

APPENDIX: USEFUL BOUNDS 

Lemma A.l. Let, for i in some finite index set I, ah, be sequences such 
that 

(19) limsup — loga^*-* < Wj. 



Then 



n— >oo n 



limsup — logfy^al*-*^ < := max 
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Proof. Although we beheve the proof is rather standard, we present it 
here. Choose an arbitrary e > 0. Then (19) entails that there is an rii such 
that for all n> rii we have that On'* < exp(n(a;j Hence, for n > maxj Tij, 

a'^'^ <exp(n(a;* + e)). 

Then 

limsup-logf^a^,^') < limsup - log({#/} x e"('^*+^)) = w'' + e. 

The stated follows after sending e | 0. □ 

Lemma A. 2. (i) There exist positive constants k and A such that 

P(M >x)< Kexpi-Xx^-^"). 
(ii) There exists positive constants R and A such that 

oo 

exp(-fc2-2^^) < Acexp(-Air2-2H)_ 

k=K 

Proof. Part (i) follows immediately from Duffield and O'Connell [12], 
Section 3.2; part (ii) is due to Dieker and Mandjes [11], Lemmas 2.1 and 
2.2. □ 
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